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ABSTRACT 

The uniformity of the helium-to-hydrogen abundance ratio in X-ray emitting intracluster medium 
(ICM) is one of the commonly adopted assumptions in X-ray analyses of galaxy clusters and cos- 
mological constraints derived from these measurements. In this work, we investigate the effect of 
He sedimentation on X-ray measurements of galaxy clusters in order to assess this assumption and 
associated systematic uncertainties. By solving a set of flow equations for a H-He plasma, we show 
that the helium-to-hydrogen mass ratio is significantly enhanced in the inner regions of clusters. The 
effect of He sedimentation, if not accounted for, introduces systematic biases in observable properties 
of clusters derived using X-ray observations. We show that these biases also introduce an apparent 
evolution in the observed gas mass fractions of X-ray luminous, dynamically relaxed clusters and 
hence biases in observational constraints on the dark energy equation of state parameter, w, derived 
from the cluster distance-redshift relation. The Hubble parameter derived from the combination of 
X-ray and Sunyaev-Zel'dovich effect (SZE) measurements is affected by the He sedimentation process 
as well. Future measurements aiming to constrain w or Hq to better than 10% may need to take 
into account the effect of He sedimentation. We propose that the evolution of gas mass fraction in 
the inner regions of clusters should provide unique observational diagnostics of the He sedimentation 
process. 

Subject headings: diffusion — cosmological parameters — X-rays: clusters 



1. INTRODUCTION 

Clusters of galaxies are powerful cosmological probes 
and have the potential to constrain properties of dark 
energy and dark matter. Recent development in X-ray 
observations of galaxy clusters have produced a large 
statistical sample of clusters a nd start to delive r pow- 
erful cosmological constraints (lAllen et al.l 12004 120081 : 
iMantz et all 120081 : IVikhlinin et alJ l2008f > that are com- 
plimentary to and competitive with other techniques 
(e.g., supernova, baryon acoustic oscillation, and weak 
lensing). This has motivated construction of the next- 
generation of X-ray satellite missions (e.g., eROSITA) 
to push the precision cosmological measurements based 
on large X-ray cluster surveys. However, in the era of 
precision cosmology, the use of clusters as sensitive cos- 
mological probes require solid understanding of cluster 
gas physics, testing of simplifying assumptions, and as- 
sessing associated systematic uncertainties. 

One of the commonly adopted assumptions in X-ray 
cluster analyses include the uniformity of the hclium- 
to-hydrogen abundance ratio with nearly primordial 
composition in X-ray emitting intracluster medium 
(ICM). At present, there is no observational test of 
this assumption, since both H and He in the ICM are 
fully ionized, which makes it difficult to measure their 
abundances using traditional spectroscopic techniques. 
Theoretically, on the other hand, it has long been 
suggested that heavier He nuclei slowly settle in the 
potential well of galaxy clusters and cause a concen- 
tration of He toward their center (|Abramopoulos et al.l 
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Ettori fc Fabianl l2006fh In the era of precision cosmol- 
ogy, this could be a source of significant systematic 
uncertainties in X-ray measurements of galaxy clus- 
ters and cosm ological paramet er derived from these 
measurements (|Markevitchll2007ft . 

Thus, the primary goal of the present work is to as- 
sess the validity of this assumption and associated sys- 
tematic uncertainties in X-ray measurements of key clus- 
ter properties as well as cosmological parameters derived 
from these observations. In this work, we investigate 
the effects of He sedimentation on X-ray measurements 
of galaxy clusters by solving a set of diffusion equations 
for a H-He plasma in the ICM. By taking into account 
observed temper ature profiles obtained by recent X- 
ray observations (IVikhlinin et al . 2006; P ratt et al.ll2007l : 
iLeccardi fc MoleridlT2008l : iGeorge et al.l 120081) . we show 
that the observed temperature drop in the cluster out- 
skirts lead to a significant suppression of He sedimenta- 
tion, compare d to the results based on the isothermal 
cluster model ([Chuzhov fc Loebll2004h . Our analysis in- 
dicates that the He sedimentation has negligible effect 
on X-ray measurements in the outer regions of clusters 
(e.g., rsoo), and it does not affect cluster mass mea- 
surements obtained at the sufficiently large cluster ra- 
dius. The effect of He sedimentation, on the other hand, 
introduces increasingly larger biases in X-ray measure- 
ments in the inner regions and could affect cosmological 
constraints, including the dark energy equation of state 
parameter w derived from distance-redshift relation as 
well as Hq derived from the combination of X-ray and 
Sunyaev-Zel'dovich effect (SZE). 

The paper is organized as follows. In § [21 we describe 
the dependence of X-ray clusters measurements on He 
abundance in the ICM. The physics of He sedimentation 
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and cluster models are discussed in §[3] In §|4l we present 
results of our He sedimentation calculations and investi- 
gate their effects on cluster properties and cosmologi- 
cal constraints derived from X-ray cluster observations. 
Main conclusions are summarized in § [6] 

2. EFFECT OF HE ABUNDANCE ON X-RAY 
MEASUREMENTS 

The observed X-ray surface brightness of a cluster is 
primarily from bremsstrahlung continuum emission of 
electrons scattering off of protons and He nuclei, which 
is given by the integral of the emission along the line of 
sight, 



S x oc dl (n e n p A ep + n e n H eA e He) 
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where x = nn D /n p is the He-to-H abundance ratio and 
A ei is the band-limited cooling function resulted from 
free- free emission of electrons scattering off ion species i, 
which is proportional to Zf. The derivation of the sec- 
ond expression takes into account the charge dependence 
of cooling function and the fact that the number density 
of electron is given by n e — n p + 2nn c — n p (l + 2x). 
Since Sx is observed and fixed, the proton number den- 
sity n p inferred from X-ray observations depends on x, 
as n p oc 1/^/(1 + 2x)(l + Ax). The derived gas mass 
density therefore depends on x as, 
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The gas mass M gas is the volume integral of equation ([3|) . 
The hydrostatic mass profile of a spherically-symmetric 
cluster depends on the local value of He abundance 
through the mean molecular weight of particles /i, 
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Note that the hydrostatic mass at a fixed mean over- 
density, A, is given by Ma oc (T e //i) 3 ' 2 . The gas mass 
fraction is then defined as / gas = M gas /M to t- 

For the primordial abundance (X = 0.75, Y = 0.25 by 
mass), the He-to-H abundance ratio is x = 0.083 and 
the mean molecular weight is /i = 0.59. If, for example, 
the Hc-to-H abundance ratio is enhanced by a factor of 
two from the primordial value, the analysis based on the 
assumption of the primordial abundance causes an un- 
derestimate of /? gas by 5% and an overestimate of M to t 
by 12%. 

3. HELIUM SEDIMENTATION IN X-RAY CUSTERS 

3.1. Diffusion Equations 

Particle diffusion in clusters is characterize d by the 
Burg ers equations of a multicomponent fluid (|Burgersl 
1969). Each species s obeys an equation of continuity 
and momentum conservation, 
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Here the species s has mass A s m pi charge Z s e, density 
n Sl partial pressure P s , and velocity u Sl where m p is 
the proton mass. The center-of-mass of a fluid element 
moves with a velocity u = ^ s n s A s u s / n s A s . The 
differential, or diffusion, velocity between species s and 
the fluid element is then w s = u s — u. These diffusion 
velocities satisfy mass and charge conservation, 



A s n s w s = 

s 

Z s n s w s =0 
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Note that the summations include both ions and elec- 
trons. To satisfy these conservation laws, for each sink- 
ing helium nuclei, there are roughly four protons and two 
electrons that float up. 

Equation ([6|) describes forces acting on a species s, 
and it is the balance of these forces that ultimately de- 
termines the rate of sedimentation. For a sinking He nu- 
cleus, the gravitational force (<?) is counteracted by three 
types of forces provided by the induced electric field (E), 
the pressure gradient (dP s /dr) of the ICM, and the drag 
force due to collisions with surrounding particles. Note 
that we neglected small terms including the inertial term 
(du s /di) and the shear stresses (or viscosity) due to col- 
lisions among the same species. We also neglected terms 
related to the coupling of thermal and particle diffusions, 
which lead to an underestimate of the He-to-H mass ra- 
tio by 20% in the radial range considered in this work 
(see §Ej|. 

Note that the sedimentation destroys hydrostatic equi- 
librium since redistribution of particles introduce a tem- 
poral change in the total gas pressure. However, hydro- 
static equilibrium can be restored quickly. This equilib- 
rium restoring acquires a net inflow with a mean velocity 
u, 
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where P gas = n s k&T is the total gas pressure for ideal 
gas, and pg as = ^ s n s A s m u is the gas density. Equa- 
tions ©-© describe the process of particle diffusion in 
the ICM. 

3.2. Resistance Coefficient 

When the plasma is sufficiently rarefied, which is the 
case for the ICM, particle pairs interact via a pure 
Coulomb potential. In the absence of magnetic field 
and turbulence, the resist ance coefficient is given by 
IChaTi man fc Cowlind {1952) as. 
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where /i st = A s A t m p / (A s + A t ) is the reduced mass of 
species s and t, and a typical value of the Coulomb log- 
arithm is lnA st ~ 40 for the H-He plasma in the ICM. 
Resistance coefficient describes the momentum transfer 
rate of species s per unit volume due to collisions with 
species t (in units of gcm~ 3 s _1 ). The resistance coeffi- 
cient, K st oc T -3 / 2 , is inversely proportional to the ICM 
temperature. Therefore, particle transport is more effi- 
cient in the ICM with higher temperature. Note that 
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heat transport via the same particle collision physics 
gives that the thermal conductivity depends on T as 

kocT 5 / 2 . 

However, magnetic field and turbulence, present in real 
clusters, can significantly modify the resistance coeffi- 
cients. To date, theoretical work on this subject has pre- 
dicted a wide range of the magnetic suppression factor, 
Jb = k/ks p ~ 0.1 — 1, wh ere ks p is the Spitzer thermal 
conductivity (|Spitzeij|1962f ). depending sensitively on the 
strength and the geometry of magnetic field as well as 
the nature of MHD turbulence. In the case of a tangled 
magnetic field, the thermal conductivity may be moder- 
ately s uppressed (fs ~ 0-1 — 0.2) relative to the Spitzer 
value (jNaravan fc Medvedevl 120011 : IChandran fe Maronl 
I2004T) - Magnetothermal instability tends to drive mag- 
netic field lines to a ra dial direction, which gives the sup- 
pression of !b > 0.4 dParrish et al.ll2008h . It was also 
suggested that MHD turbulence may provide the mag- 
netic suppression fac tor of order unity (|Cho fc Lazarianl 
[2^0llLa^a7ia^l2W7h . 

Observationally, one can use the strength of observed 
temperature gradients in the ICM to constrain the ef- 
ficiency of thermal conductivity. Observed large-scale, 
negative ICM temperature gradients set the upper limit 
on the thermal conductivity, fs < 0. 2 for a Tx = 10 keV 
cluster with a typical age of 7 Gyr (|Loebl 12003 ). Obser- 
vations of A754 also yield a suppression fac tor fs < 0.1 
for the bulk of the ICM (|Markevitchll2003h . The width 
of cold fronts observed by recent X-ray observations sug- 
gests a_ronsiderablv_jm^er sup_pression factor (fs <C 
0.1) (jMarkevitch fc Vikhlinm1l2007l for a review). This, 
however, is a local constraint which is not applicable to 
a cluster as a whole. One might also imagine that the 
abundance profile of heavy nuclei (e.g., Fe) with X-ray 
emission lines may provide further insights into the effi- 
ciency of particle diffusion. Though possible in principle, 
interpretation of the abundance profile is complicated by 
the fact that the ICM is enriched continuously by strip- 
ping of metal-enriched gas from cluster galaxies, and that 
the process of sedimentation is generally slower for heav- 
ier nuclei. 

Given the current large uncertainties in diffusion coeffi- 
cients due to the lack of knowledge on magnetic field and 
turbulence in the ICM, we parametrize the effective resis- 
tance coefficients as K st = f^Kf^ , where we take Jb 
as a free parameter. In order to illustrate how our results 
depend on the suppression factor, we consider two cases 
with Jb — 1 (un-magnetized) and 0.2 (tangled magnetic 
field). The former should be taken as an extreme model, 
while the latter is roughly corresponds to the current ob- 
servational limit discussed above. Note, however, that 
all current observations are consistent with fg — (con- 
duction fully suppressed). We, therefore, caution that 
fs may be orders of magnitude below the unity. 

3.3. Sedimentation Velocity 

To develop physical insights into the process of He sed- 
imentation, it is useful to consider a drift velocity of 
a trace He particle in a background of hydrogen, i.e, 
n p ^> nHo and w p — 0. In this limiting case, equa- 
tions ([6|) for the two species are decoupled. The right- 
hand side of the equations for H vanishes, thereby fixing 
an electric field, eE — 0.5 m p g. Substituting E into the 
equation of motion for He, we obtain the sedimentation 



velocity of He nuclei as WHe = 3 m p griYi e / K p yi c , which 
gives 

iuhc — 80 kms _1 /s 

J—Y /2 ( I \ ( n p N-i 

)keVy VlO- 7 - 5 cms- 2 MlO- 3 cm- 3 / V J 

where the induced electric field counteracts gravity and 
suppresses the sedimentation speed by 25%. At a fixed 
density, the sedimcntion velocity is generally larger for 
higher temperature (T) and gravity (g). Note that the 
pressure gradient (dP/dr) term in equation (|6|) further 
suppresses the sedimentation process. Typical sedimen- 
tation timescale in clusters is generally longer than the 
Hubble time. The equilibrium distributions to the Burg- 
ers equation (lAbramopoulos et al.lll98lHQin fc Wull2000l : 
iChuzhov fc Nusseril2003r)~ are therefore not applicable for 
clusters. Instead, a full time-dependent calculation is re- 
quired for this analysis. 

3.4. Cluster Models and Initial Conditions 

We set up cluster models and initial conditions as fol- 
lows. Initially, we assume that the ICM consists of a 
primordial H and He plasma uniformly throughout clus- 
ters. We ignore the contribution of heavier elements of 
Z > 2. 

For the total mass distribution, we adop t the Navarro- 
Frenk - White (NFW) density profile (|Navarro et al.1 
[1993), 



p{r) 
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x'{\ 



-A 2 



(12) 



where x' = r/r s , r s is a scale radius, and p s is a nor- 
malization constant. Mass enclosed within a radius r is 
then given by M(x') = 47rp s r 3 [ln(l + x 1 ) - x'/(l + x')]. 
Throughout this work, we define the cluster mass to be 
Ma = (4-7r/3)r A A p cr n, where ta is a radius of a spheri- 
cal region within which the mean enclosed mass density 
is A times the critical density of the universe p cr it- We 
adopt A = 500 and the concentration C500 = ?"500 /r s = 4 
(V06). We also consider A = 2500, where some of the X- 
ray measurements are also made. In this cluster model, 
»"250o/?*500 = 0.46. 

We consider two ICM temperature models: (1) the 
isothermal temperature profile T(r) = Tx, and (2) the 
observed temperature profile obtained with deep Chan- 
dra observations of nearby relaxed clusters (V06) given 

by 

T{r) _ 1 (i/0.045) 1 9 + 0.45 1 
~Tx~~ ' (x/0M5y- 9 + l [1 + (5/0.6) 2 ] ' 45 ' 

(13) 

where x — r/r^Q, and Tx is the X-ray spectral tem- 
perature, which is related to M 50 o = M 5 (Tx/5 keV) Q , 
where M 5 = 2.89 x 10 14 h- 1 M Q and a = 1.58 (V06). 
For an illustration, we plot these ICM temperature pro- 
files in Figure[TJ The observed temperature peaks around 
r /f500 — 0.2, and decreases at both inner and outer radii. 
For example, the temperature drops by nearly a factor 
of 2 from the peak value at r/r^ = 0.01 and 1. 

We set up the initial gas distribution by assuming hy- 
drostatic equilibrium of the ICM in the potential well of 
clusters dominated by dark matter, 

kB dM ^ -P^9(r) , (14) 
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Fig. 1. — Temperature profiles of X-ray clusters used in our 
models: isothermal temperatur e profile (dashed line) a nd observed 
temperature profile obtained by[Vikhlinin ct al. (2006) (solid line). 

where the gravitational acceleration at a radius r is 
g(f) = GM(r)/r 2 . For each temperature profile T(r), 
we derive the initial gas density profile as, 



Pgas(»") — Pgas,0 



To 
T(r) 



exp 



9(r') 
T(r') 



dr' (1,5) 



where the central density p gas ,o is normalized by requir- 
ing M gas /M(r 5 oo) = /gas, 500 = 0.15 (V06), and T is the 
central temperature for each ICM temperature model. 
For isothermal case, the density profile reproduces the 
analytical form given by p sas (r) = jo gas ,o exp(-rjfj,)(l + 

x^M x ' with r\ = ATTGm p p s r 2 J{k B Tx) (|Makino et all 
fl998h . 

Finally, we consider two types of mass accretion histo- 
ries (MAHs) of clusters: (1) a static cluster with a fixed 
mass, (2) the MAHs of cluster-size halos calibrated with 
N-body simulations of cluster formation. For the latter, 
we adopt an a nalytical expression o f averaged MAHs of 
halos given by Ivan den Boschl (|2002l ). 



-0.301 



log(l 



log(l + z f ) 



(16) 



where Mq is the cluster mass at present time, Zf is the 
formation redshift defined as M(zf)/M — 0.5, and v 
is a parameter that determines the shape of MAHs of 
dark matter halos, and it is strongly correlated with Zf. 
This functional form provides a reasonable description of 
the MAHs of simulated c lusters for a typical va lue of v 
in the range of 1.4 — 2.3 tTasit siomi et al.ll2004f) and an 
average formation redshift of (zf) ~ 0.6 ( Cohn fc White 
120051 see also Berrier et al. 2008). 4 The points that 
are particularly relevant for this work are that a typical 
present-day cluster grows by a factor of two in mass in 

4 Note that the equation J 16H provides a reas onable descrip- 
tion of the MAHs of simulated galactic-size halos l|Wechsler et al.l 



2002) by tuning a free parameter / to 0.254 (see Appendix A 
in Ivan den B osch 2002). For the cluster-size halos, we find that 
/ = 0.656 provides a reasonable description of their MAHs and the 
average formation time of a large statistic al samples of clusters- size 
halos extracted from N-body simulations (Cohn & White 2005, sec 
also Berrier et al. 2008). 



Fig. 2. — The spatial distribution of helium-to- hydrogen mass 
fraction (Y/X) in a Tx = 10 keV static cluster with fg = 1. Lines 
indicate two types of ICM temperature profiles shown in Fig. [T] 
From bottom to the top, the set of curves shown are for cluster 
ages of 0, 1, 3, 5, 7, 9, and 11 Gyr, respectively. 

the past « 6 Gyr since z f=a 0.6, while an equal-mass 
cluster at higher redshift (z = 1) forms in a consider- 
ably shorter timescale (w 1-2 Gyr). This suggests that 
the high-z clusters are generally much less affected by 
the process of He sedimentation than their low-rcdshift 
counterparts. In §4.1-4.3, we use the simple static cluster 
model to investigate the dependence of the He sedimen- 
tation efficiency on some of the key cluster parameters, 
including the temperature, the age, and the magnetic 
suppression factor. In § 14. 4[ we use the realistic MAHs 
to investigate evolution in X-ray observable properties 
with redshift. Throughout this work, we use cosmologi- 
cal para meters: 17m = 0-3, = 0.7, — 0.0462, and 
h = 0.7 (|Komatsu et al.ll2008ft . 

4. RESULTS 

4.1. Spatial Distribution and Evolution of He 
Abundance 

Starting from the initial cluster model described in 
§ 13.41 we follow the diffusion process in the H and He 
plasma by solving equations ©-© numerically. At each 
time step, equations ©-© are solved to obtain diffusion 
velocities (w s ) for H, He and electrons as well as electric 
field (E). Using these diffusion velocities, we update the 
abundance of each species by solving equations (O and 
([9]) together. We repeat this procedure through the mass 
accretion histories of galaxy clusters. We use 600 spatial 
grids logarithmically spaced in a computational domain 
of 10~ 3 < r/r 50 o < 10 and set the diffusion velocities to 
zero at both inner and outer boundaries. These choices 
ensure that results are robust in regions of our interest 
(0.01 < r/r 500 < 1). 

Figure [2] shows evolution of a spatial distribution of 
the He-to-H mass fraction (Y/X) of a T x = 10 keV 
cluster with Jb = 1. Here we compare results of the 
isothermal model and those of the V06 model. This 
comparison shows that the observed temperature drop 
in the outer regions significantly suppresses the process 
of He sedimentation in the outskirts (r/r^vi > 0.3) of 
clusters. The efficiency of He sedimentation, however, is 
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Fig. 3. — Biases in X-ray measurements of galaxy clusters at 
^2500 as a function of Tx for a static cluster age of 7Gyr. Lines 
indicate the isothermal model (dashed line) and the V06 model 
with fg = 1 (thick-solid line) and fg = 0.2 (thin-solid line). 

enhanced in the inner regions near the peak of the ob- 
served temperature profile, and it is suppressed again in 
the inner most regions (r/rsoo < 0.02) due to the ob- 
served temperature drop at these radii. At r = rsoo, 
the He abundance in the V06 model is very close to the 
primordial value (Y/X = 0.333). The effect of He sedi- 
mentation is negligible at this radius, and these results 
are fairly independent of cluster age and the value of /b ■ 
These conclusions are in stark contrast with the result 
of the isothermal model which gives considerably larger 
He abundance of Y/X = 0.479. In the V06 model, the 
value of Y/X increases rapidly toward inner regions and 
peaks at r = 0.06 rsorj. The mass fraction ratio increases 
to Y/X = 0.405 at r = r25oo and 1.0 at r = 0.1 rsoo for 
the cluster age of 7 Gyr. 

4.2. Effect on Gas Mass, Total Mass, and Gas Mass 
Fraction 

As we discussed in § [3J the spatial variation in the 
He abundance caused by the mass segregation could in- 
troduce observational biases in X-ray measurements of 
galaxy clusters. Figure [3] shows biases in X-ray measure- 
ments of gas mass, total mass, and gas mass fraction 
at r25oo for three sedimentation models: the isothermal 
model with f B = 1 (dashed line), the V06 model with 
/b = 1 (thick- solid line) and /s = 5 (thin- solid line). Re- 
sults are shown for a typical cluster age of 7 Gyr and for a 
range of average cluster temperatures (Tx = 1 — 10 keV). 
Here we define these biases to be ratios between the quan- 
tities derived by assuming the primordial abundance and 
their true values in our cluster models, which corresponds 
to the values that would have been obtained if the sed- 
imentation effect has been taken into account in X-ray 
data analyses. 

Comparing results of two different models, we find that 
biases in X-ray measurements of total cluster mass at 
^2500 are considerably smaller in the V06 model. In the 



Fig. 4. — Biases in X-ray measurements of galaxy clusters r25oo 
as a function of cluster age. The curves starting from the flat line 
corresponds to cluster ages of 0, 1, 3, 5, 7, 9, and 11 Gyr. Results 
are shown for static clusters of the V06 model with fg = l. 

case of a hot (Tx — 10 keV), un-magnetized (/b = 1) 
cluster, for example, the total mass is overestimated by 
4% and 14% for the V06 model and isothermal model, 
respectively. The biases in M gas , on the other hand, are 
underestimated by 4% for both models. Note that these 
values arc similar for two different temperature models. 
This is because Af gas measurements, obtained by inte- 
grating the gas density over the cluster volume, depend 
primarily on the average ICM temperature (Tx), and 
relatively insensitive to the details of ICM temperature 
profiles. / gas is then biased low by 8% and 16% for V06 
model and isothermal model, respectively. Note that the 
bias in / gas above is obtained by integrating from the 
cluster center to the radius r. This yields a large bias 
than the one evaluated by using the local He abundance 
at r (|Markevitchll2007T ). 

The biases in X-ray measurements are very sensitive to 
the average ICM temperature (Tx) as well as the mag- 
netic suppression factor (/b). The biases are consider- 
ably smaller in cluster with smaller values of Tx or /b, 
because the sedimentation process is slow in these clus- 
ters with larger resistance coefficient. We find that the 
biases at r 2 5oo become negligible (< 2%) for Tx < 3 keV 
or /b < 0.2. These biases also depend sensitively on 
the age of clusters. Figure Q] illustrates that these biases 
increase rapidly with the cluster age and cause the appar- 
ent evolution of X-ray measurements of hot, dynamically 
relaxed clusters. Our sedimentation model based on the 
observed V06 temperature profile predicts that the bi- 
ases in / gas are likely less than 10% for a realistic range 
of cluster parameters (including the cluster age, Tx, and 
/b). Note also that the biases become negligible in clus- 
ter outskirts (r ~ r§oo) for all systems. 

4.3. X-ray and SZE-derived Hubble Constant 

Measurements of angular diameter distance derived 
from the combination of X-ray and Sunyaev-Zel'dovich 
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Fig. 5. — Evolution of cluster gas mass fraction enclosed within 
^2500 as a function of redshift. Res ults are for a 10 keV cluster with 
realistic MAHs given by eq. 1161 . Lines indicate the isothermal 
model with fg = 1 (dashed line), V06 model with fg = 1 (thick- 
solid line) and fg = 0.2 (thin-solid line). Dotted lines indicate the 
evolution in / gas caused by changes in the dark energy equation of 
state w by 10% relative to the fiducial ACDM cosmological model 
with no sedimentation (indicated with a straight dotted line). 



Effect (SZE) observations also depend on the assumed 
He abundance, x, as, 



d,A oc 



y 2 l + Ax 
S X T? 1 + 2x 



(17) 



where the expression was obtained by canceling elec- 
tron densities in both the X-ray surface brightness Sx oc 
n 2 e d A (l + 4x)/(l -I- 2x) (see eq. [2]) and the SZE comp- 
tonization parameter y oc n e T e d,A- Since Hq oc d^ , 
the X-ray+SZE derived H measurements could be af- 
fected by the He sedimentation process. As shown in 
the bottom-right panels of Figures [3] and 2J our model 
suggests that the sedimentation could introduce biases 
in the X-ray+SZE derived Hq measurements high by 
about < 15%, with the exact value depending on the 
cluster age, temperature, as well as the magnetic sup- 
pression factor. This effect is of the same sign and or- 
der as the 6% offset seen between the X-ray and SZE- 
derived Hubble constan t (H = 77.6±|;| kms" 1 Mpc" 1 ) 
(|Bonamente et al.l [2006) and res ult of the Hubble Ke y 
Project (H = TSkms" 1 Mpc" 1 ) (jFreedman et alJ2001f ). 
Given the errors in the current X-ray+SZE Ho measure- 
ments, our models with fg = 1 and 0.2 are both consis- 
tent w ith the observed offset. As proposed bv lMarkevitchl 
(|2007h . the comparison of the X-ray+SZE derived Hq 
and independent Hq measurements could be used to con- 
strain the He abundance in clusters. 

4.4. Evolution of Cluster Gas Mass Fraction 

The evolution of cluster gas mass fraction of X-ray lu- 
minous, dynamically, relaxed clusters can provide pow- 
erful observational constraints on the equation of state 
of dark energy w 



Sasaki 1996; Pe n! Il997t lEttori et al 
2003|lAllen et al.ll2004tlLaRoaue et alJl2006tlAllen et al 
20081) . The sensitivity to w lies in the dependence of the 



Fig. 6. — Evolution of cluster gas mass fraction enclosed within 
three cluster-centric radii, r = 0.1 rsoo (solid lines), r = 0.3 rsoo 
(dashed lines), and r = r25oo (dotted lines). R esult s are for a 10 
keV cluster with realistic MAHs given by eq. Q2), for the V06 
model and fg=X. 

observed / gas to the angular diameter distance, cIa(w, z), 

d A (w,z) 



f gas (w, z) = 0.113 



d\(w = —1.0, z) 



1.5 



(18) 



where the observed / gas is 0.113 at z = 0. Recent 
measurements of the / gas evolution based on Chandra 
X-ray observations of 42 bright, dynamically relaxed 
galaxy clusters yielded w = —1.14 ± 0.31, by assum- 
ing a flat geometry and standard priors on flbh 2 and 
h. The combined analysis of / gas plus CMB and SNIa 
measurement s constrains w to better than 10%, w = 
-0.98 ± 0.07 (|Allen et al.l[2008h . 

Here we point out that the effect of He sedimentation 
can introduce an apparent evolution in X-ray measure- 
ments of cluster gas mass fractions, which could lead to 
systematic biases in observational constraints on the dark 
energy equation of state, w. This is illustrated in Fig- 
ure [SJ where we show evolution of the observed gas mass 
fraction at r25oo by using the realistic, time-dependent 
MAHs given by equation (fl6|) in order to follow the 
He sedimentation process in a growing cluster-size halo, 
starting from a proto-cluster at z = 5. Our sedimen- 
tation models show that the observed gas mass fraction 
is overestimated by 12%, 6%, and 2% at z = 1, for the 
isothermal model with /b = 1, V06 model with fg = 1 
and with Jb = 0.2, respectively. The apparent evolution 
m /gas arises because clusters at lower-redshifts have had 
more time to experience the sedimentation on average. 

As shown in Figure [3 the effect of He sedimentation is 
degenerate with changes in the dark energy equation of 
state, w. Note that the changes in w by -10% (relative to 
the fiducial ACDM cosmology with w = — 1) corresponds 
to the changes in / gas by about 3% between z = and 1. 
The effect of He sedimentation is to make the best-fit w 
more negative, i.e., w < — 1. For example, results of the 
V06 model with /g = 1 (0.2) are degenerate with the cos- 
mological models with w = —1.18 (—1.04) without sed- 
imentation. The isothermal model with fg = 1, which 
should be taken as an extreme case, requires w = —1.38. 
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Given the measurement uncertainties (~ 30% in w) , cur- 
rent constraints on w should not be affected significantly 
by the effect of He sedimentation. However, future X-ray 
measurements aiming to constraint w to better than 10% 
will need to take into account of this effect. Note that 
the effect would be larger for a population of dynamically 
relaxed clusters than the average population considered 
here. 

Given that the variation in the He abundance in clus- 
ters are unknown at present, it would be interesting to 
ask whether the current X-ray data could be used to 
constrain the abundance distribution. Here we propose 
that the / gas evolution in the inner cluster regions, where 
the effect of sedimentation is expected to be larger, can 
be used as sensitive probes of the He sedimentation pro- 
cess. In Figure HI we illustrate evolution of / gas (V06 
model with = 1) at three different radii, including 
r = O.lrsoo, 0.3r5oo, and r25oo- This figure shows that 
changes in / gas between z = and 1 are amplified signif- 
icantly in the inner regions of clusters. At r = O.lrsoo, 
/gas evolves by ~ 30% from z = to 1, which could be 
detected with current data. If the gas mass fraction evo- 
lution is detected, it can provide unique constraints on 
the efficiency of He sedimentation, which in turn provides 
constraints on the magnetic and turbulent suppression of 
particle diffusion in the weakly magnetized ICM. Non- 
detection is also interesting, as it will provide an upper 
limit on the efficiency of particle diffusion in the ICM. 

5. DISCUSSIONS 

Here we comment on additional effects that may af- 
fect our sedimentation calculations. First, turbulent 
mixing will likely play an important role in determin- 
ing the efficiency of sedimentation. Turbulence, on the 
one hand, will tend to mix fluid elements from differ- 
ent parts of the cluster and to counteract the effect of 
sedimentation, but it may also enhance the rate of sed- 
imentation by increasing the mobility of the ions. De- 
tailed investigation of the effect of turbulence on sed- 
imentation process is out of the scope of the present 
work. But, we point out that recent hydrodynamical 
simulations of galaxy cluster formation uniformly indi- 
cate the pre sence of ubiquitous subsonic turbulent flow 
in clusters (Norman & Brvan 1999; Nagai et al. 2003; 
Rasia et ail |2004 iKav et alJ 12004 iFaltenbacher et al 
20051 iDolag et alj 120051 : iRasia et all 120061 : iNagai et al 
2003). These simulations indicate that the turbulence 
can provide about 5 — 10% of the total pressure sup- 
port at r = rsoo even in relaxed systems, and its relative 
importance increases in cluster outskirts as well as unre- 
laxed clusters with recent major mergers (Lau, Kravtsov, 
Nagai, in preparation). It is therefore possible that tur- 
bulent mixing caused by gas accretion in cluster outskirts 
and/or major mergers may significantly modify the ef- 
ficiency of the sedimentation. Hydrodynamical cluster 
simulations that also include the sedimentation physics 
will likely provide realistic assessment of this effect. 

We also point out that the standard diffusion equa- 
tions assume that the mean-free-path of the ions is small 
compared to the size of the system. The proton mcan- 
free-path, however, could become a sizable fr action of 
the cluster size in their outskirts (jLoebl [20071 ). For a 
Tx = 5 keV isothermal cluster, for example, the mean- 
free-path of proton is comparable to the cluster virial 




500 



Fig. 7. — The ratio of two He sedimentation velocities derived 
with (lojje) an d without (i«He o) the effect of thermal diffusion, 
plotted as a function of the cluster-centric radius. 



radius (w 2 x r^oo). The validity of the fluid approxi- 
mation, however, depends sensitively on the ICM tem- 
perature profile. For the cluster with the same Tx, the 
mean-free-path in the V06 model, for example, is consid- 
erably smaller (9%) fraction of the virial radius, which 
makes the fluid approximation still reasonable in the ra- 
dial range of our interests. 

In this work, we also neglected the effect of thermal dif- 
fusion, which can cause the change in the diffusion speed. 
To assess the effect of thermal diffusi on, we follow the 
formalisms presented in Burgers (1969) which takes into 
account the heat flow interferences, and apply it to the 
H-He plasma in the ICM. For isothermal case, the diffu- 
sion velocity is enha nced by ~ 23%, which a grees with 
the value reported bv lChuzhov fc Loe b (2004). For non- 
isothermal case, we find that thermal diffusion enhances 
He sedimentation in the outskirts where the temperature 
gradient is negative, but suppresses He sedimentation in 
the inner region (r/r^oo < 0.1) with positive temperature 
gradient. Our analysis indicates that the contribution 
of thermal diffusion to the diffusion velocity by about 
< 50% throughout clusters. Figure [7] shows the change 
in He sedimentation speed caused by the effect of thermal 
diffusion as a function of r/Vsoo- Note that the thermal 
diffusion causes relatively small change in the Y/X at 
cluster outskirts (e.g., r ~ r 50 o), despite its large frac- 
tional change in diffusion velocity. It is because that the 
diffusion speed was very small there to start with. Using 
the sc aling of Y/X with time given by IChuzhov fc Loebl 
(|200l . we find that the effect of thermal diffusion on 
Y/X is less than 20%. 

It has been suggested that the enhancement of He 
abundance in the central regions might cause the declin e 
of the observed Fe abundance (Ett ori fc Fabianl l2006f) . 
To assess this effect, we repeat the He sedimentation 
calculations for the three cool-core cl usters (includin, 
Centaurus, A2199, A1795) presented in lEttori fc Fabia 
(2006), using their observed gas density and temper- 
ature profiles. Our calculations indicate that the en- 
hancement of He in the core of the Centaurus cluster is 
less than 20% from the primordial He abundance, even 



Peng & Nagai 



in the extreme case of un-magnetized cluster with the 
cluster age of 11 Gyr. This corresponds to a decline 
of the observed Fe abun dance by only 7% (see Fig. 4 
of lEttori fe Fabiar]|2006ft . Similar results are found for 
A2199 and A1795. Note that the He sedimentation is sig- 
nificantly suppressed in the cool core regions. We thus 
conclude that the observed 20 — 50% reduction in the 
iron abundance in some of the cool-core clusters cannot 
be explained by the He sedimentation alone. 

6. CONCLUSIONS 

In this work, we investigate effects of He sedimentation 
on X-ray measurements of galaxy clusters and their im- 
plication for cosmological constraints derived from these 
observations. By solving a set of flow equations for a H- 
He plasma and using observationally motivated cluster 
models, we show that the efficiency of He sedimentation 
is significantly suppressed in the cluster outskirts due to 
the observed temperature drop, while it is dramatically 
enhanced in the cluster core regions. Our sedimentation 
model based on the observed temperature profile sug- 
gest that the effect of helium sedimentation is negligible 
a t r 500j and it does not affect cluster mass measurements 
obtained at the sufficiently large cluster radius. However, 
the effect of sedimentation increases toward the inner re- 
gions of clusters and introduces biases in X-ray measure- 
ments of galaxy clusters. For example, at T2500 7 biases 
in X-ray measurements of gas mass, total mass, and gas 
mass fractions, are at the level of 5 — 10%. The effect 
of He sedimentation could also introduce biases in the 
estimate of the Hubble parameter derived from the com- 
bination of X-ray and SZE measurements, which could 
explain the observed offset in the X-ray+SZE derived 



Hq and independent measurement from the Hubble Key 
project. We emphasize, however, that the magnitude of 
these biases depends sensitively on the cluster age, tem- 
perature, and magnetic and/or turbulent suppression in 
the ICM. 

We show that the process of He sedimentation intro- 
duces the apparent evolution in the observed gas mass 
fractions of X-ray luminous, dynamically relaxed clus- 
ters. The effect of He sedimentation could lead to biases 
in observational constraints of dark energy equation of 
state w at a level of <10%. These biases tend to make 
the value of w more ne gative. Current measu rements 
based on / gas evolution \ Allen et all 12004 12008ft should 
not be significantly affected by these biases. However, fu- 
ture measurements aiming to constrain w to better than 
10% may need to take into account the effect of He sed- 
imentation. For cosmological measurements, one way to 
minimize these biases is to extend the X-ray measure- 
ments to a radius well beyond r25oo- At the same time, 
the evolution of cluster gas mass fraction in the inner 
regions of clusters should provide unique observational 
diagnostics of the He sedimentation process in clusters. 
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